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1 Introduction 



Strange particle production in heavy ion collisions is a subject of great interest since it may 
provide evidences of the formation of a new state of matter: the quark-gluon plasma (QGP). 
The tantalizing pursuit of this discovery is the justification for the efforts being made to carry 
out the heavy-ions collision experiments that have taken place at the AGS and CERN-SPS, 
and to plan new ones for future facilities. A phase transition between normal nuclear matter 
to a deconfined phase has been predicted long ag dii to occur at very high temperatures and 
densities. For the time being, the properties of the QGP have not being clearly stablised, 
neither the order of the phase transition between the two states. This makes the problem of 
recognizing the formation of QGP, if it happens, particularly difficult, more so since we have to 
rely on indirect clues to its detection. Several signatures have being proposed for this purpose 
that may deliver circumstantial evidences of the formation of quark matter (for an up-to-date 
review on the subject with a complete list of references, see 0). The strange mesons, baryons 
and antibaryons production has been predicted to be greatly enhanced in the presence of QGP 
as compared to the case when there are only standard hadron matteiS i . Also the study of the 
stopping power may help to uncover the presence of the deconfined phase: a rapid change in 
the shape of the rapidity baryon distributions with increasing incident beam energy is thought 
to be a clear signal for new degrees of freedonii. 

An additional problem is caused by the lack of a full description of the confined phase, which 
constitutes the background to the events in which QGP may have been formed. Enhanced 
strange particle production is a meaningless expression unless we can find out what strange 
particle production must be expected from an ordinary ion- ion collision event. We can not 
extract this information from QCD since we do not know how to apply this theory when the 
effective coupling constant is not small. This is also the ultimate cause of the lacunae in 
our understanding of the QGP properties and of the order of the phase transition mentioned 
in the above paragraph. We have to rely on phenomenological models to fill this gap. The 
Dual Parton Model (DPM)! and the Quark Gluon String Model (QGSM)i i have been very 
successful in describing many features of hadron-hadron, hadron-nucleus and nucleus-nucleus 
collisions. Nevertheless, they were not able to reproduc^M the rapidity distributions of A and 
A hyperons, as well as proton minus antiprotons {p — p), in central sulphur-sulphur collisions 
measured by the NA35 Collaboration0' @ . The idea that the onset of the deconfinement 
phase transition may be at the origin of this disagreement is very appealing. However, we must 
carefully explore any other possible explanation of these data before admitting that eventuality. 
The NA35 data show a large production of A's in the fragmentation as well as in the central 
region. It is very difficult to understand why the formation of QGP in a collision between 
identical nuclei should contribute so largely to the borders of the accessible phase space. This 
can not obviously be due to QGP formation as it should mainly affect the hadron spectra in 
the central rapidity region. Besides, QGP formation in this experiments, involving relatively 
light nuclei, would come as a striking surprise. The energy densities currently available at 
the AGS and CERN-SPS of 0.5-10 GeV/fm^ correspond to temperatures between 100 and 200 
MeVEl, for which the effective coupling constant of QCD is not less than 1, and we should not 
be expecting a deconfinement phase transition to occur. This also should encourage us to look 
for conventional explanations to the phenomena observed in heavy nucleus collisions that have 
not been understood so far. 
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Strangeness enhancement is also present in pp interactions. The ratio K/ tt increases both 
with the energy of the colhsion and with the multiphcity^l . It has also been observed an increase 
of the ratio A/p in pp collisions as the centre-of-mass energy increases from 20 GeV to 1800 
GeVli. Again, the NA35 Collaboration also reported that strangeness is enhanced between pp 
and pS collision^l. At the same time, the average number of collisions per participant nucleon 
increases between pp and pA collisions as well as between pA and central AB collisions. This 
increase is also correlated with an increase of the multiplicity of secondaries. According to the 
DPM and the QGSM, the extra particles are produced by the fragmentation of strings that 
have sea quarks at their ends. As this mechanism on its own failed to reproduce the observed 
enhancement of A and A in central SS collisions, additional sources of this increase have been 
explored. 

Diquark-antidiquark pairs from the nucleon sea have been introducecf0'0lil. This model 
predicts, contrary to the experimental evidence, the same increase in the absolute value of 
strange baryons and antibaryons. It may also be worth mentioning that this mechanism is 
very similar in its formulation and its results to the, so called, string fusion model of ref. . 



In reference |]TU| the problem is treated considering final state interactions of co- moving pions 
and nucleons of the type tc + N K + A. This approach is not convincing since it contains 
the unrealistic assumption that all baryons are created instantaneously at the beginning of the 
collision so the secondary collisions can take place. 

What the data on strange baryon production of the NA35 Collaboration may actually be 
suggesting is that we may be missing an important piece in our understanding of the nucleon- 
nucleon rescattering. I devote this paper to reconsider how the valence and sea quarks take part 
in the formation of the ends of strings. I will show that, by allowing the leading diquarks to be 
formed not only of valence quarks, a good understanding of the data in the fragmentation region 
is obtained. I also considered the net proton distribution, usually regarded as a measurement 
of the stopping power of nuclear matter. Data form the NA35 Collaboration shows a larger 
stopping than predicted. The model proposed here also gives better results for this particular 
problem. 

There have been a couple of other mechanisms proposed that could contribute to the full 
explanation of the data. One of thenS assumes that diquarks can break up into two quarks, 
one of which can be slowed down along the stringcarrying the string junction (baryon number) 
and produces much slower baryons. The otheiH consists in considering the possibility that 
the baryon number is not always carried by a diquark, or a quark as in the previous case, 
but by a non-perturbative configuration of gluon fields. Both have similar phenomenological 
implications as both mainly affect the centre of the rapidity space and can not account for what 
happens in the fragmentation region what is my major concern in this paper, in particular for 
the case of A production. 

There are a number of other models that are also constructed around the idea of independent 
fragmentation of strings. The best known is the Lund-Fritiof modelil . Although the general 
assumptions are at first sight very similar to the ones in DPM and QGSM, the theoretical basis 
from which these models were developed are quite different and I will not refer to them in the 
rest of the paper. 

In section ^ I discuss the theoretical basis of the paper and make a full description of the 
model. The next two sections deal with fits to pp data to obtain the value of some parameters 
of the model, and the results for nucleus-nucleus collision. The paper ends with the conclusions. 
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Figure 1: Representation of a cut reggeon and its relation with the planar diagram. 

2 The Model 

The independent-fragmentation string models, DPM and OGSMjWere inspired by the devel- 
opments of the Dual Topological Unit arizat ion (DTU)© eI 0' 0. The amplitudes of high- 
energy hadronic collisions were formally written as an expansion in powers of l/Nf, where Nf 
is the number of flavours. Each term in the expansion can be associated with a diagram with a 
definite topology. The scheme was based on the S-matrix analysis and was closely related to the 
Regge Field Theory (RFT). The DTU program was reinforced when Veneziancfil established 
a conceptual link with QCD as he included in the series powers of l/iV,§'0i), being the 
number of colours. In this way, the whole series of the topological expansion formally includes 
all the terms of the perturbative QCD series ensuring unitarity. 

The lowest order term (fig. |l| ) is the planar diagram connecting two colliding-hadron lines. 
It represent the exchange of an object that has the same quantum number of the Reggeon and 
can be identified with it. Unitarity requires the inclusion of all diagrams obtained by repeatedly 
multiplying the planar term by itself. The next order topology we have to consider is that of 
the cylinder that corresponds to the Pomeron (see fig. |^) as it only vacuum quantum numbers 
can be exchange in the elastic collision. Higher order terms can be obtained by combining the 
simplest topologies described above, and all of them have a one to one correspondence with 
Reggeon Field Theory (RFT) diagrams. Bearing in mind that we are interested in the high- 
energy regime, we make use of this correspondence to neglect terms containing the exchange 
of one or more reggeons: although they may be dominant in the topological expansion, their 
contribution decreases with the centre-of-mass energy squared s as s°'^^^^~^ ~ s~^'^ while the 
pomeron exchange subdiagrams behave like s'^^'C')^! = s^, where A > 0. 

This formal developments laid the foundation of phenomenological models that had to fill 
the conceptual gaps left by the theory. Both DPM and QGSM postulated that multiparticle 
production in hadronic collisions at ISR energies, take place in two definite steps. First, the 
hadrons split in two fragments that share the total momentum of the hadron according to the 
corresponding structure functions. One of the fragments belongs to the 3 representation of the 
SU{3) colour group, a quark, and the other to the 3 representation, an antiquark in the case of 
a meson or a diquark in the case of a baryon. Secondly, two colour strings are formed each one 
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Figure 2: Diagrammatic representation of a cut pomeron and its relation with the cyhnder 
diagram. 

having a fragment of the 3 representation in one end from one of the hadrons, and another of 
the 3 representation, from the other hadron at the other end. These two strings decay into two 
chains of hadrons that form the final multiparticle state. The contribution of a string to the 
final state is given by the fragmentation functions that are, along with the structure functions, 
the other basic ingredient of the models. 

When we move to higher energies, we have to introduce terms with more than two strings 
that correspond to diagrams with several cut pomerons. The probability weights of the different 
diagrams are given by the perturbative RFT that can be effectively described by a generalized 
eikonal modefl . Each cut pomeron add two chains of hadrons to the final state as the result 
of the decay of two more strings. The generalization of the models to hadron-nucleus and 
nucleus-nucleus collisions carries some similarities with this extrapolation to higher energies. 
In this case, multiple inelastic collision appears as a consequence of the rescattering of each 
nucleon by the nucleons on the other interacting nucleus. Again, for each additional inelastic 
collision we have to consider the fragmentation of two extra chains. The probability weights 
are given in this case by an external model like that of Glauber-Gribov, that reproduces the 
geometry of the collision. 

Both DPM and QGSM assume that the diquarks at the end of the strings in hadronic 
collisions, are always made of valence quarks. So, in the case of a colliding proton, regardless of 
the number of cut pomerons in the collision, the diquark can only be of type {uu) or (ud). The 
central point to this paper is the assumption that the sea quarks can be part of the diquark 
with the same probability as the valence quarks. It is still true that the diquark in a proton 
can only be of the types mentioned above when it suffers only one inelastic collision. But, when 
the number of cut pomerons is greater than one, we also have to consider diquarks of the type 
{dd), (us), (ds) and (ss). The flavour content of the diquark depends on the flavour content of 
the hadron sea. 

This assumption means the inclusion of terms of the topological expansion that are missing 
in the other independent strings fragmentation models. They correspond to diagrams where 
the cylinder (the pomeron) is linked to the colliding hadron through a q'"q'^ or q"q'^ trajectory 
as depicted in fig. 0(b). For the incoming nucleon to interact with two nucleons of the target 
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Figure 3: Two diagrams corresponding to inelastic collisions of a nucleon with two nucleons 
from the target. Diagram (a) is common to all models, while diagram (b) is one of the new 
components considered in this paper. 

it is necessary that it appears in a virtual state composed of a baryon-like and a meson-like 
states. Then, the topology of the interaction can be depicted as the superposition of a baryon- 
nucleon and a meson-nucleon collisions. In other models, the baryon-like state could only be 
formed with the three valence quarks, and only meson-like states of the type \uu >, \dd >, 
etc, where considered. I am dropping this restriction only justified for high-mass intermediate 
states, where sea quark- ant iquark pairs are created with large values of rapidity (> 1) respect to 
their parent hadron. Therefore, in a proton - nucleus collision we now have to include, besides 
\p* 7r°*(p, ■■■)>, states like |A* K+* >0. 

These new components enhance the production of the strange baryons in the fragmentation 
region due to the contribution of strange diquark fragmentation. I consider a theory with only 
three flavours: u, d and s, and assume that the u and d contents of the sea are equal, while the 
s content is suppressed and represented by the parameter S. 

- . a) 



u + d 



SU(3) symmetry corresponds to a value = 1. I take S = 0.5 that is compatible with several 
data analysidi'i'ii. In the rest of the paper, I will be using 7 as the portion of strange 
sea-quarks over the total number of sea-quarks: 



7= ' . (2) 

u + d + s S + 2 



For the given value of S, 7 = 0.2. 



^ Baryon-like and meson-like, as well as the asterisks, mean that they are not real baryons or meson since 
only the flavour content is considered and no attention is paid to other quantum numbers 
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2.1 Structure Functions 



The only essential difference between the DPM and the QGSM comes from the way they 
postulate the structure functions for sea quarks and antiquarks. Only the leading behaviour 
at a; —> is retained in all cases, that has been called the retardation of the hadron fragments. 
In a collision in which there is only one cut pomeron, the retardation of a quark of flavour / 
depends on the intercept of the corresponding q-^q^ trajectory while for the diquark depends 
on the intercept of the corresponding qqqq exotic trajectory. When dealing with nucleons, for 
these functions we have 

Uix) = f,,{x) = . (3) 

Where we have made use of the equation, vahd in the planar approximation, aqaqb^aqb + aqcqc = 
'2oiqaqbqc and ^^(O) (aAr(O)) is the zero intercept of the leading Regge (nucleon) trajectory. 

When there are more than one cut pomeron, we need quarks and antiquarks from the sea to 
make up enough string ends. In the DPM, it is assumed that quark-antiquark pairs are formed 
as the result of the decay of a pomeron which has a retardation given by 1/x. Following this 
line of reasoning DPM postulates that the retardation of a sea quark or sea antiquark is given 
by 

Where /x is a parameter that sets a mass scale. QGSM postulates, on the other hand, that the 
retardations of the sea quarks and antiquarks have the same origin as that for valence quarks 
and, hence, are the same functions. I adopt this latest view in the present work. While there are 
no practical difference in terms of computational results at ISR energies, I find QGSM structure 
functions more coherent with the assumption I have made - sea-valence diquark symmetry. The 
mechanism considered in DPM to postulate the retardation functions, must be included in a 
treatment of the problem in which we sum over all possible RFT diagrams. For instance, in the 
case of two inelastic collisions, one term must correspond to a coupling between the hadron and 
the two cut pomerons mediated by the exchange of a pomeron (DPM) - PPP diagram, another 
one to the same coupling mediated by the exchange of a reggeon (QGSM)- PRP diagram. The 
corresponding relative weights for this kind of diagrams are not known and their computation is 
an open problem. The choice among neglecting one or the other constitutes a model assumption 
that can only be justified by the phenomenologyS . Yet, the eikonal approximation to RFT 
points that the PPP diagrams constitute the main contribution to high-mass intermediate 
states and eventually lead to the point-like quarks structure functions, while the PRP type 
of diagrams are the dominant terms for low- mass intermediate stateS. This idea allowed to 
understand the data from HERA on small-x behaviour of the structure functiondSl. 

The PPP coupling does not allow the inclusion of diagrams of type 0(b). The extra sea 
quark-antiquark pairs, in collisions with more than one cut pomeron, are the result of the decay 
of an object with vacuum quantum numbers. Therefore, DPM only allows diagrams of type 
0(a). This also explains why only valence diquarks contribute to the final state in DIS and 
hard hadron-hadron collisions. Then, as far as we are interested in the low-j9^ regime, it seems 
most reasonable to adhere to QGSM view. 

The s-quark, (us), (ds) and (ss) retardations are given by the following functions, where 

is the intercept of the ip trajectory and Aa = a_R(0) — a^(0). 
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fs{x) 

fus (3^) 
fssix) 




.Qfl(0)-2aAr(0)+2Aa 



,afl,(0)-2ajv(0)+Aa 



(5) 



The corresponding functions for the uu, ud, and dd diquarks, as well as for the u and d quarks 



The momentum-fraction distribution function in a nucleon, in the case of cut pomerons, 
is given by the following expression: 



This function is perfectly defined if we consider that fi{x) = fi{x). The indexes ij = {u, s, d} 
and ij = {u, d, s}, id = {{uu), (ud), (dd), (us), (ds), (ss)}. Since we are considering nucleons as 
the colliding particles, the diquark can be of type {ss) only if > 3 and of types (us) or {ds) 
if iV > 2; we find the s-quark retardation only when N > 2. 

The single-parton structure functions are obtained integrating over all other momentum 
fractions: 



are given in equation ^. 




(6) 




C'o''x--^{l-xf{l-^xf-\ N>1 



P. 






X [1 -7 + 7(1 
Co^^''^x(l - xf-kl - -fxf~\ N > 1 



(7) 




(jqq,N^aR{0)-2aN(0)+Aa^^ _ ^ya^{0)+2(N -l){l-aR{0)) 

X [1 -7 + 7(1 -xf^T"' 
Cr'^xi(l-x)^-i(l-7x)^-2, N>2 
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The numerical values for the constants are ai?(0) = 0.5, Aa = 0.5 and = —0.25. The 
value of ajsf has been chosen within the theoretically allowed range to obtain the best agreement 
with the proton-proton data. 

The factors [1 — 7 + 7(1 — x)^^"]^""' in equation (|^) have their origin in the fact that the 
s quarks take a larger portion of the nucleon momentum that the u or d quarks. It has a little 
effect on the structure functions when N is small. When is large, these factors make the 
functions steeper at a; — > 0. 



2.2 Fragmentation Functions 

The fragmentation functions are derived from Regge argumentdS' S . The ones needed for the 



present work are listed in the appendix. I have corrected some missprints in ref. I^H] and use 
consistently the value of given aboveQ. The constant A = 2a'p']_ ~ 0.5, where a' is the slope 
or the Regge trajectory, has been used in the expressions. 

All these functions are proportional to one of the constants and a^, which repre- 

sent the universal rapidity densities of protons, A hyperons, antiprotons and anti-A hyperons 
respectively, around the centre-of-mass of the string. The fragmentation of a diquark into a 
given baryon contains two terms, one corresponds to the creation of a hadron that contains the 
baryon number of the diquark: primary or leading fragmentation (subindex 1), the other comes 
from baryon-antibaryon formation in the string: secondary fragmentation (subindex 2). Hence, 
the latest is proportional to the corresponding antibaryon density. Functions such as Di y^^(z), 
include a factor function linear with z. This is introduced to take into account the excess in the 
production of baryons containing the whole diquark over those that contains only one of the 
quarks of the diquark. Taking into account the constants that define these factors, cq, ci and C2, 
we have seven parameters with only three of them independent as a consequence of the baryon 
number sum rule. I have assumed that this sum rule saturates with only protons, neutros and 
A's as if they were the only stable baryons under strong interactions, being all the strong decays 
taken into account by the Regge treatment of the problem. In fact there are other hyperons 
that decay through electroweak processes into A's and nucleons which have some contribution 
to the final rapidity spectra. We can easily include in the treatment these hyperons assuming 
that the corresponding fragmentation functions are proportional to the ones for the A produc- 
tion. This is strictly true for the cases of S"*",S° and S~ while it is only approximate for the 
case of S hyperons which provide a much smaller relative contribution. I take the contribution 
from this decays making use of the available phenomenological information. The total number 
of hyperons (anti-hyperons) is taken to be 1.6 the number of A's (A). 30% of the other strange 
hypeons (anti-hyperons) decay into A (A)0' § . Hence, with all these considerations taken, we 
may interpret a\ (a^) as the hyperon (anti-hyperon) density in the middle of the string some of 
which will decay into A (A) and others into protons (antiprotons) and neutrons (antineutrons) 
in equal proportions. The corrections introduced by these consideration are not significant, in 
the case of A distributions, given the accuracy of the experimental data as will be seen latter 
in the following sections. 



■^The treatment of the fragmentation functions in ref. |10[ is inconsistent since two different values for a^, 
—0.25 and —0.5, were used simultaneously. This second value was also used to compute the structure functions. 

^In ref. [|To[ the constraint imposed by baryon number conservation was overlooked and different values for 
the baryon densities in the strings were used. 
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The functional form of the diquark fragmentation functions for diquarks with strangeness 



-1 and —2 were obtained applying the results of ref. [^. I have assumed that both strange and 



non-strange quarks are dynamically symmetric, hence the factor [1 + ^/z) /2. The fragmentation 
of a diquark into a baryon that do not contain any quaxk from the diquark is taken to be 0. 
This may be controversial since it has been suggested^ that an excess in the production of 
Vt~ over fi"*", in non-strange hadrons collisions, contradicts this assumption. As a matter of 
fact, preliminary data show a very small asymmetry in vr-nucleus interactions at 500GeV/cii'. 
Nevertheless, although the data are still inconclusive, this effect, if small as it seems, may be 
understood in the framework of the model described in this paper: multiple collision raises the 
probability of having contributions to the final state coming from the fragmentation of strange 
diquarks that results in a net production of VL~ . 



2.3 Rapidity Distributions 

We concentrate in the study of collisions between identical nuclei. The atomic weight and the 
atomic number are A and Z respectively. The rapidity distribution for the production of a 
hadron h is given by: 



dy (Tab 



{ns - UAm^^A^f (y) + ^A^fiy)) + - ^i^)(<;f (?/) + (?/))] (§) 

-|-Symmetric terms (n^ ^ ns)} 

In writing this expression, I neglect kinematical correlations among different strings, 
Here, ^""^'"S'^'A'A'/^s is the cross-section for AB collisions where there are n inelastic nucleon- 
nucleon collisions involving ua nucleons of A and ns nucleons of 5, in which at least one of 
the participating nucleons of A (B) suffers fiA if^s) inelastic collisions. 

I restrict the discussion to the case A = B. We can approximate equation ||, with a good 
degree of accuracy using average values for ua, and n. 



dN 



= ^mi--^^ {y) + {y)) + (n - ns){Nl^Z iv) + ^1^1 (^)) • 0) 

Equation is obtained by neglecting the dependence of the string density on the /x's, 
which are replaced by their average values p, 

KulM = E^WNJJ^dx^dx2g%ix,)g%ix2)^—iy-A,sr''^) (10) 



e{sXiX2 - si' 

Here a and b refer to quark or diquark. A sum is performed over all possible partons flavours 
i and j. The probability weights and are trivially computed from the model assumptions. 

^In practical computations a slightly more accurate expression has been used in which the average of densities 
over = were taken using the weights o^j^j a a- These weights are obtained from ref. No significant 
numerical difference were observed compared to the results obtained using /I^ ~ njuA- 
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They depend on 7, A and Z. A mass threshold for each chain is enforced by the ^-function. 
Please note that this mass threshold is not always equal to the mass of hadron h since the 
flavour structure of the string many times forces to have other hadrons besides h in the final 
state . A is the center-of-mass rapidity of the string. 



;i2) 



dy "'^'^ '~\Gi^{y) \iy<A 

where G{y) are the fragmentation functions: 

Gl{y) =Z^D^^iZ^) 

Here Z+ = exp[(y - A) - yuAx] and Z_ = exp[-(?/ - A) - yMAx], Vmax being the maximum 
rapidity that hadron h can have in the string for fixed values of momentum fractions xi and 

X2. 



3 Proton-proton collisions 

Since I have not considered sum rules others than the ones due to baryon conservation, the 
model leaves three free parameters that I have fixed in order to obtain the best fits to proton- 
proton data. At the ISR energies the structure of these collisions is very simple as we only need 
to consider one pomeron exchangd0. This makes the choice of this energy range particularly 
convenient as we avoid having to take into account the probability weights for the different 
re-scattering configurations within the same nucleon-nucleon collision. The only differences 
between the results of the present model and those from DPM and QGSM, at this stage, come 
from the choice of a^, that was usually taken in the past to be near —0.5. 

Figure ^ represents the Feynman x distribution of protons produced in the reaction pp 



pX. The fit obtained agrees very well with the experimental data taken from ref. partic- 
ularly in the region xp < 0.75. It has been pointed out that experimental points above that 
cut are not very reliable (see [^) due to detector inefficiencies. The fit in the region x ^ is 
very sensitive to the value of a^. For a^r = —0.5 the results of the fit are much poorer that 
the one shown, which corresponds to the value of = —0.25 consistently used throughout 
this paper. The value of Op fixes the values of the remaining parameters by the use of the sum 
rules. Figure ^ shows the fit to pp —>■ pX. It is also worth mentioning that the same value for 
ap was independently obtained from both fits. Please note that the computation performed do 
not include diffraction and, consequently, the results obtained do not show the characteristic 
diffraction pick at large-x. 



To obtain de value of a^, I used the data from ref. that were taken at a similar centre- 
of-mass energy. The reactions considered are pp AX and pp AX. We show the best fits 
in figures || and |^ (please mind that in this two fits only is free, since oa is given by the 
sum rules and the value of ap) The large error bars in the experimental data do not allow for 
a very accurate determination of a^. Corrections from contributions of electroweak decays of 
heavier hyperons, as discussed in subsection p.2| , have also been considered. I have taken the 
value that best fits the data and also that is compatible with the mass suppression factor of 



introduced in [^| 
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The numerical values obtained are: Op = 0.667, qa = 0.282, ap = 0.077, = 0.057, 
Co = 1.430, ci = 6.076 and C2 = 67.560. 

4 Nucleus-Nucleus collisions 

At this point we have a model for nucleus-nucleus collisions which contains no free parameters. 

The result for the rapidity distribution of A hyperons produced in central sulfur-sulfur 
collisions is compared with the experimental data from the CERN-NA3d0 Collaboration in 
figure ^ The short-dashed line corresponds to the results of the model without the sea-valence 
diquark symmetry introduced in this paper. As we can see, both models give similar results 
at ycAi = 0, while there is a remarkable difference in the fragmentation region where only the 
present model reproduces the data. The data points are not corrected for particles originated 
from electromagnetic S° decays and from weak S (S) decay. These decays have been taken 
into accout as described in pl2[ The uncorrected results are showed with a long-dashed line 
and represent the total strange hyperon distribution. In computing the correction an effective 
recoil, due to the decay of the resonances of 0.2 units in rapidity was introduced which has none 
or very little effect on the results. Figure ^ shows the agreement between the A distribution 
and the data points from the same experiment. Again, the dashed line corresponds to the 
uncorrected hyperon distribution. The only disagreement is with the point near i/cm = in 
the A distribution. 

Figure shows the proton minus antiproton distribution in central (solid line) and pe- 
ripheral (dashed-dotted lirie) sulfur-sulfur collisions. The data points were also taken by the 
CERN-NA35 collaborationEJ. For central collisions, the agreement in the fragmentation region 
of the spectrum is very good while the data show an overpopulation of protons around ycM = 0. 
The data were extracted from the measured charge asymmetry of the final state applying model 
dependent corrections. The cancellations oi p — p and 7r+ — tt^ have been assumed due to the 
isospin symmetry of the system. The results showed here do not make use of the first assump- 
tion above. A small difference, not significant in magnitude, is obtained in the central rapidity 
region when we compute the leading proton distribution as if the p — p cancellation was exact, 
leaving the integrated correction equal to zero as it should. From quark flavour counting the 
model predicts an overal 7r+ — tt^ cancellation and an negligible change on this result due to 
any local effect. Corrections based on a Monte Carlo simulation (Fritiof) to the data have 
also been applied to estimate the — K~ contribution to the charge asymmetry. The model 
presented here predicts a much higher production of A, in agreement with the data, than Fritiof 
whose results are very close to those of this model without the new components (dashed line 
fig. P). For each extra produced lambda we are left with a s quark at the end of one string that 
will appear in the final state as a or a in equal proportions. The obvious conclusion 
is that this model predicts a higher — K~ asymmetry concentrated in the central rapidity 
region and whose total value is equal to one half the excess of A predicted by this model over 
the predictions of Fritiof. I have added a Gaussian correction to the p — p distribution that 
accounts for the missed point in the central region (dashed line). Given the long error bars, 
the precise shape of this correction do not matter significantly. It is interesting to notice that 
another correction to the A distribution might be necessary due to the effect on the A — > pvr" 
background coming from — *■ tt+tt". Although, there is not enough information available to 
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Figure 6: Fit to the Feynman x differential cross section da/dx (given in mb) for the inclusive 
production of A in the reaction pp — > AX. Data correspond to p^^^^ = 405 GeV and were taken 
from ref . ^ . 
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Figure 7: Fit to the Feynman x differential cross section da/dx (given in mb) for the inclusive 
production of A in the reaction pp AX. Data correspond to Pi^}^ = 405 GeV and were taken 
from ref . |^ . 
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Figure 8: Rapidity distribution of A's produced in central sulphur-sulphur collisions represented 
in the nucleon-nucleon centre-of -mass reference frame. The solid line correspond to the results 
of this model while the short-dashed line are the results obtained when the new components 
are not included. The long dashed line represent the total strange-hyperons spectrum. The 
data points are from ref. ||TT| and have been reflected onto positive rapidities. 
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Figure 9: Rapidity distribution of A's produced in central sulphur-sulphur collisions represented 
in the nucleon-nucleon centre-of-mass reference frame (solid line) and the total strange hyperons 
spectrum (dashed line) . The data points are from ref. [|lT| and have been reflected onto positive 
rapidities. 
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assume that this background has not been correctly subtracted in the experimental analysis. 

Figure ^ shows the predictions for the rapidity spectra of A and A, solid lines, produced 
in central Pb-Pb collisions. The experimental points correspond to preliminary data from 
the NA49 CoUaborationil. The dashed lines correspond to the total hyperon (antihyperon) 
distributions. Figure |12| compares the result for the net proton distribution with data from the 
same collaborationEl . These data points were obtained subtracting the background due to the 
decays A pvr" using the results of a Monte Carlo simulation. Correcting for the difference 
between the A spectrum predicted by this model and the one used in the experimental analysis, 
we obtain the solid line out of the net distribution (dashed-dotted line). I have also corrected 
for the — K~ asymmetry as explained above, dashed line. As we can see, the experimental 
data tend to show a larger stopping in Pb-Pb collisions that what the model predicts. 

As mentioned before, I have made use of the Galuber-Gribov model to obtain the average 
numbers of participating nucleons and nucleon-nucleon collisions. I have used the value cx]^ = 
30 mb for the inelastic nucleon-nucleon cross section. The profile functions where obtained using 
the nuclear matter distribution given in ref. |4^. They consist of three parameters, individually 
fitted to nuclear data, functions and are more accurate than the Saxon- Wood density functions 
used in the past. The actual numerical values used for central S-S, = 27.8 and n = 62.8, 
were obtained requiring that b < 1, where b is the impact parameter. For peripheral collisions 
we have ua = 9.7 and n = 14.7, while for central Pb-Pb collisions are ua = 204.4 and n = 895.9. 

To perform the computations, I have taken the values of the proton and A transverse mass 



equal to 1.1 and 1.3 GeV respectively, in agreement with the data from [15 



5 Conclusions 

In order to explain the data on baryon rapidity spectra presented by the NA35 collaboration, 
I have included into the DPM and the QGSM models components that were previously over- 
looked. They are compatible with the 1/A^ expansion of QCD on which these two models are 
based. A discussion over the right structure functions for the string's ends leads to the adoption 
of the ones proposed by the QGSM, as they are the only ones consistent with the model's as- 
sumptions. The fragmentation functions are obtained from the planar approximation to Regge 
diagrams. They have been published elsewhere except for the ones corresponding to strange 
diquark that are published here for the first time. All of them are given in the appendix. Once 
taken into account the baryon-number sum-rules, the number of free parameters of the model 
is reduced to 3. The numerical values of these parameters have been obtained fitting data of 
inclusive baryon and antibaryon production in proton-proton collisions. The sum rules allow 
us to obtain the fragmentation functions that do not contribute to proton-proton collisions at 
ISR energies. 

The model has strong theoretical grounds and all the approximations done are well justified, 
being the uncertainties introduced by them small and well under control. 

At the level of nucleus-nucleus collisions, we have a model which contains no free parameters. 
The A distribution for central sulphur-sulphur collisions was computed and a perfect agreement 
with the experimental data was obtained. This result does not depend on the new components 
introduced in this paper; all the A particles come from baryon-antibaryon formation in the 
strings, which is essentially independent of the parton components at the end of the string. This 
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Figure 10: p — p rapidity distribution in sulphur- sulphur collisions represented in the nucleon- 
nucleon centre-of-mass reference frame. The dashed-dotted line corresponds to primary protons 
produced in peripheral collisions; the solid line are the results for the proton minus antiproton 
distribution in central collisions. The data points are from ref. [0 and have been reflected onto 
positive rapidities. The dashed line represents the correction for the excess of — K—. See 
text for a more detailed explanation. 
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Figure 11: Rapidity distribution of A's and A's, solid lines, produced in central PbPb collisions 
represented in the nucleon-nucleon centre-of-mass reference frame. The dashed lines represent 
the total hyperons (anti-hyperons) distributions. Data points are from ref. 
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Figure 12: Rapidity p — p distribution in central PhPh collisions represented in the nucleon- 
nucleon centre-of-mass reference frame. The solid line corresponds to the prediction including 
the correction from hyperon electroweak decays into protons applied to the initial proton dis- 
tribution, dashed-dotted line. The dashed line represents the correction for the excess K'^ over 
K~ (see text for an explanation). Data points are from ref. E3. 
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result is very relevant since it clearly shows that the s content at the central rapidity region is 
very well understood with no need to assume collective effects such as QGP formation. 

I have also computed the rapidity distribution for the A hyperons in central sulphur-sulphur 
collisions and compared it with data from the same collaboration. The results show a remarkable 
increase of the rapidity distribution in the fragmentation region well in agreement with the 
experimental data. This is due to the inclusion of chains initiated by diquarks that can contain 
a quark from the nucleon sea. The fact that this sea-quark has a probability of being strange, 
obtained by independent analysis, is responsible for this effect. The only disagreement with 
the data is seen at the central rapidity region. 

The net proton distribution for the same collisions has also been obtained and compared with 
data of the same collaboration. We observe a good agreement between the model results and 
the experimental points in the fragmentation region. Once the correction of the electric charge 
asymmetry due to the excess of over K~ is included, we also understand the spectrum at 
mid-rapidities. 

The NA35 data on baryon distributions do not support that the formation of QGP has taken 
place or that we need to take into account collective effects of any sort. I have compare the 
model results with the preliminary data on A and A rapidity distributions in Pb-Pb collisions 
form the NA49 Collaboration; a partial agreement with the data was found. Unfortunately there 
are no points available in the fragmentation region. The comparison with the p — p spectrum 
measured by the same collaboration shows that the model under-predicts the stopping power 
in Pb-Pb collisions. This raises the interesting question about which of the model assumptions 
fails for the case of very heavy nuclei. It could be that effects like those of references and 
[pT| have a significant contribution to this stopping, falling within the experimental errors in the 
case of sulphur-sulphur collisions. This possibilities should be careful check in hadron-hadron 
interactions as they contain the implicit assumption that the wave function of the colliding 
hadron couples to the pomeron over a rapidity region larger than one unit. This may have 
an effect on diffractive scattering. The assumption that collisions between heavy nucleus can 
be decomposed as nearly uncorrelated nucleon-nucleon collisions, as implied by the use of the 
Galuber-Gribov model, may also fail when the size of the nuclei involved becomes very large. 

Another important conclusion, that has been obtained as a byproduct, is that the analy- 
sis performed on the proton-proton data strongly support a value for the intercept of the N 
trajectory close to —0.25. This is also a major result of the model that in return can help to 
establish a more accurate value for the mass of some resonances. 

For future works it is left the detailed study of the meson sector. Also the implications of 
the new structure of the pomeron-hadron coupling discussed in the paper will be investigated 
in order to shed light on other outstanding problems such as heavy quark formation in hadronic 
and nuclear collisions. 
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Appendix 

I list here the fragmentation functions into baryons and antibaryons. The ones not shown can 
be obtained applying isospin symmetry p n, u d. 
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